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ABSTRACT 


A finite element analysis of the end notch flexure specimen for mode II 
interlaminar fracture toughness measurement was conducted. The effect of 
friction between the crack faces and large deflection on the evaluation of 
Gjic from this specimen were investigated. This analysis showed that the 
friction effect is negligible when crack length is greater than one-fourth 
of length of the specimen. Further, geometric nonlinear analysis is required 
to evaluate Gj I( , when deflection is greater than the thickness of specimen. 
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INTRODUCTION 


The delamination of laminated fiber reinforced composites is of 
current interest, particularly with regard to their durability and 
damage tolerance. It is the most commonly observed failure mode 
in composites. Fracture mechanics has been found to be a useful 
tool for understanding delamination [1]. Most of the work reported 
in this area has been concerned with the delamination in the opening 
mode, i.e. mode I. Delamination under shear sliding mode i.e. mode II 
and combined mode I — 1 1 situations are also the possible fracture modes 
of potentially greater significance in composite structures. Several 
studies on mixed-mode interlaminar fracture are also reported [1-4], 
Recently attention has been given’ to mode II shear delamination 
[5-9], Several investigators have used end notch flexure (ENF) 
specimen to measure shear mode strain energy release rate, G||£, for 
delamination growth [6-9]. This specimen is a beam with a crack 
(delamination) located on the neutral plane at one end and is 
subjected to the three-point bending. This specimen was initially 
utilized by Barrett and Foschi to measure Mode II critical-stress- 
intensity factor, K||£, for a strongly anisotropic fibrous material, 
wood [10]. Currently an ASTM task group is conducting a round 
robin test program using an end notch flexure specimen as a prelude 
to developing an ASTM standard for measurement of mode II interlaminar 
fracture toughness of composites. In this study, the end notch flexure 
specimen was analyzed using the finite element method for understanding 
its behavior and performance. 
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FINITE ELEMENT ANALYSIS 


The end notch flexure specimen geometry is shown in Fig. 1. 

Two geometries of the specimen were analyzed; (1) total length, 

2L = 7.62 cm and (2) 2L = 10.16 cm. For both geometries, the 
thickness (2h) of the specimen was taken equal to 0.30 cm which 
is nominal thickness of a 24 ply composite panel. These dimensions 
of L and h were selected since they were employed in experimental 
works [6,9]. The width of specimen was 2.54 cm. Only one material 
system was analyzed; unidirectional graphite/epoxy. The material 
properties used in finite element analysis are presented in Table 1. 

A two-dimensional finite element program called GAMNAS [11] 
was used. A typical finite element mesh, shown in Fig. 2 consisted 
of about 1000 isoparametric four-node elements and had about 2400 
degrees of freedom. The element size in the vicinity of the crack 
was 0.02 x 0.02 mm in order to evaluate accurately the strain energy 
release rate G j j . This size was selected by previous experience [12] 
as well as by conducting the convergence test on G j ^ calculation with 
mesh refinement. The error in G j j obtained in the present study is 
estimated to be less than ±2%. The strain energy release rate, Gj., 
in the analysis was computed using a virtual crack closure technique 
[13]. Plane-strain condition was assumed in the finite element analysi 

The analysis indicated that mode II deformation is achieved at 
the crack tip and that the accompanying mode I deformation causes 
closure and overlapping of the opposite faces of the crack. This is 
not possible physically. To prevent the overlapping of crack faces, 
the nodal coupling technique, available with GAMNAS program [11], 
was used. For this purpose, the multipoint constraint was applied 
at corresponding nodes to have the same displacements normal to the 



crack faces. This resulted in pure mode II condition, i.e. Gj = 0. 
When the crack faces are pressed together as in this specimen, 
friction effects will occur. The friction effects on Gjj calculation 
will be rather difficult to assess in the practical situation, however 
a reasonable estimate of its role can be investigated with the finite 
element analysis. The friction effect was modeled in the following 
manner. The analysis was conducted intially with the frictionless 
crack faces where displacements normal to the crack faces were 
matched at corresponding nodes. The normal reactions at the crack 
faces were obtained from this first step of calculation. The friction 
forces at the nodes on the crack faces were calculated for the assumed 
coefficient of friction. The finite element analysis of the specimen 
was conducted with these computed friction forces applied at the crack 
faces. This second step of analysis provided the new normal forces 
at the crack faces and corresponding friction forces. The analysis 
was repeated till the applied friction forces matched with the computed 
normal reactions. This required about two or three iterations. 
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RESULTS AND DISCUSSIONS 


Comparison Between Theoretical and FEM G | | Values 

The measured compliance, c, versus crack length, a, relationship 
of a test specimen, in general, provides the most accurate evaluation 
of fracture toughness. However, it requires several measurements. On 
the other hand, fracture toughness can be measured from a single 
measurement of critical load if the necessary expression for G calculation 
is available. The experimental dc/da values were found to be subject to 
large errors due to small variation in compliance with crack length for 
small crack length values in the ENF specimen [7J. Therefore, the 
following expression for mode II strain energy release rate, Gjj for 
this specimen was developed by Russell [7]: 



9P 2 a 2 

16b 2 E n h 3 


( 1 } 


where P, a, b, h are defined in Fig. 1 and E^ is the Young's modulus 
in longitudinal direction. The G j j , as expressed above, is based on a 
linear beam theory where shear deformation is neglected and the curvature 

1 _ d 2 y/dx 2 

R “ r , A -.2,3/2 

[1 + (dy/dx) ] 

2 2 

is approximated with 1/Rci d y/dx . The Gp values obtained from FEM 
analysis and from Eq.(l) for both geometries of the specimens; (1) 

2L = 7.62 cm and (2) 2L = 10.16 cm, are compared for the normalized 
crack length, a/L, ranging from 0.125 to 0.9. The normalized crack 
length commonly used for measurement of G j j ^ with this specimen was 
0.5 [6,7,9], Figure 3 shows that Gj|^ from a linear beam theory. 
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i.e. from Eq.(1), is within 8% and 61 of its counterpart from FEM 
analysis for a/L = 0.5 for specimens with 2L = 7.62 cm and 2L = 10.16 cm, 
respectively. Further, Fig. 3 shows that deviation between G j j from 
Eq.(l) and FEM analysis decreases with longer crack lengths and 
increases with shorter crack lengths. This comparison shows that the 
G ( | C should be evaluated from the end notch flexure specimen for 
a/L >0.5 if Eq.(l) is used. 

Effect of Overhang Length 

The effect of overhang length beyond the support was also 
investigated with finite element analysis. Specimens with overhang 
lengths of 6.35, 12.7 and 25.4 mm on each side were analyzed. The 
G|j values obtained for these cases were within ±1% of each other. 

Thus, the overhang portion of this specimen has no effect in the 
evaluation of 
Friction Effect 

As mentioned previously, the friction was also modeled in FEM 
analysis. The most difficult task is to evaluate the actual coefficient 
of friction, y, between two pressed crack faces of ENF specimen. 

This was beyond the scope of the present study. However, a simple 
experiment was performed to estimate the coefficient of friction between 
two unidirectional graphite/epoxy laminates. This experiment involved 
the measurement of the force to initiate sliding between the two 
laminates with a known compressive load. The coefficient of friction, y, 
was found to be about 0.35. Therefore, the finite element analysis 
was conducted with the three values of y = 0.1, 0.3 and 0.5 to under- 
stand the role of friction in a qualitative manner. Results of this 
analysis are shown in Fig. 4. Several observations can be made from 
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this analysis. For crack length commonly used in experimental 
investigation, i.e. a/L = 0.5, the effect of friction is negligible 
for all practical purposes if p is about 0.3. And, this effect is 
less in the ENF specimen with 2L = 10.16 cm in comparison to 
specimen with 2L = 7.62 cm. The friction effect shows a pronounced 
effect for shorter crack lengths and is negligible for longer crack 
lengths. This leads to an apparent conclusion that longer crack 
length should be employed for C|j^. evaluation, but this may cause 
the large deflection of specimen during testing and henceforth it 
may require the geometric nonlinear analysis as discussed later on. 

The friction forces between the crack faces were primarily 
in the region above the support pin. Since it is always desirable 
to have the friction force as small as possible, it might, therefore, 
be preferable to use the shims or rollers between the crack face 
just above the support pin to reduce the coefficient of friction. 

Barrett and Foschi [10] suggested the use of roller between crack 
faces for measurement of of wood from this type of specimen. 

However, it may not be convenient in composite specimens due to 
their small thicknesses. 

Geometric Nonlinear Analysis 

To investigate the effect of finite rotation due to large 
deflection, the geometric nonlinear analysis of the ENF specimen 
was conducted with the GAMNAS finite element program [11]. This 
analysis was conducted without friction. Figure 5 shows the comparison 
of Gj| obtained from linear and geometric nonlinear analysis for both geo- 
metries of specimens for crack length a = 0.5L. Linear and nonlinear solution 
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are in agreement with each other up to deflection, A , equal to the 
thickness of specimen, 2h. For larger deflection, i.e. A > 2h, the 
G|j would be overestimated from a linear analysis or Eq.(l). This 
overestimation depends on the deflection as shown in Fig. 5. 
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CONCLUDING REMARKS 


A finite element analysis of the end notch flexure specimen 
for mode II interlaminar fracture toughness measurement was 
conducted. The effect of friction between the crack faces and 
large deflection on the evaluation of G^ from this specimen were 
investigated. Two geometries of the specimen, and one material 
system was studied. This study showed that this specimen should 
be preferably used for Gj|£ measurement with crack length, 
a > 0.5L and deflection, a < 2h. 
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TABLE 1 - GRAPHITE/EPOXY MATERIAL PROPERTIES 3 



Modulus 13 , GPa 


L 

Poisson’s Ratio 0 


1 

E 22 

C 12 

v 12 

V 23 

.0 

13.0 

6.4 

0.34 

0.35 


'33 


= E 


22 ' 


13 


= v 


12' 13 


= G 


12 


The subscripts 1, 2 and 3 correspond to the longitudinal, transverse 
and thickness directions, respectively, of a unidirectional ply. 
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Fig. 3. Ratio of G j j 1 
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NORMALIZED CRACK LENGTH, a/L 

Fig, 4 (a), Ratio of Gjj with v = 0,0 to Gjj with v = 0.1 from FEM analysis 

as a function of normalized crack length. 
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Fig. 4 (b). Ratio of Gjj with v = 0.0 to G j j with v = 0.3 from FEM analysis 

as a function of normalized crack length. 
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Fig, 4 (c), Ratio of Gjj with w = 0,0 to Gjj with v = 0.5 from FEM analysis 

as. a function of normalized crack length, 
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